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Abstract 


Let G be an undirected bipartite graph with positive integer weights 
on the edges. We refine the existing decomposition theorem originally 
proposed by Kao et al., for computing maximum weight bipartite 
matching. We apply it to design an efficient version of the decomposi- 
tion algorithm to compute the weight of a maximum weight bipartite 
matching of G in O(\/|V|W'/k(|V|, W'/N))-time by employing an 
algorithm designed by Feder and Motwani as a subroutine, where |V| 
and N denote the number of nodes and the maximum edge weight 
of G, respectively and k(x, y) = log w/log(x?/y). The parameter W’ is 
smaller than the total edge weight W, essentially when the largest edge 
weight differs by more than one from the second-largest edge weight in 
the current working graph in any decomposition step of the algorithm. 
In best the case, W’ = O(|E|) where |E| is the number of edges of G 
and in the worst case, W’ = W, that is, |E| < W’ < W. In addition, 
we talk about a scaling property of the algorithm and research a better 
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bound of the parameter W’. Experimental evaluations of randomly 
generated data show that the proposed improvement is significant in 
general. 


Keywords: Weighted bipartite matching, Graph decomposition, Ex- 
perimental evaluation, Random instances of graphs, Combinatorial 
optimization. 


1 Introduction 


Throughout the paper, we use the symbols N and No to denote the sets of 
positive and non-negative integers, respectively. Moreover, the symbols N 
and W denote the largest weight of any edge and the total weight of a 
graph G, respectively. Let G = (V = Vi U va, FE, Wt) be an undirected, 
weighted bipartite graph where V; and V2 are two independent non-empty 
disjoint subsets of the vertex set V of G, and each edge in the edge set E 
connects a vertex in V; to a vertex in V2 with positive integer weight on each 
edge which is given by the weight function Wt: E > N. The (total) weight 
of the graph G is defined by W = Wt(G) = Ss Wt(e). We also assume that 


eC€E 
the graph does not have any isolated vertex. For uniformity, we treat an 


unweighted graph as a weighted graph having unit weight for all edges. 

We use the notation {u, v} for an edge e € EF’ between u € Vi and v € Va, 
and its weight is denoted by Wt(e) = Wt(u,v). We also say that e = {u,v} 
is incident on the vertices u and v; and u and v are each incident with e. 
Two vertices u,v € V of G are adjacent if there exists an edge e = {u,v} € E 
of G to which they are both incident. Two edges e1, e2 € E of G are adjacent 
if there exists a vertex v € V to which they are both incidents on. 


1.1 Fundamentals of Maximum Weight Bipartite Matching 


A subset M C E of edges is a matching if no two edges of M share a common 
vertex. A vertex uv € V is said to be covered or matched by the matching M 
if it is incident with an edge of M; otherwise, v is unmatched [2]. 

A matching M of G is called a maximum (cardinality) matching if 
there does not exist any other matching of G with greater cardinality. 
We denote such a matching by mm(G). The weight of a matching M is 
defined as Wt(M) = Wt(e). A matching M of G is a Mazimum Weight 

eC M 
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Matching (MWM), denoted as mwm(G), if Wt(M) > Wt(M’) for every 
other matching M’ of the graph G. 

Note that, if G is an unweighted graph then mwm(G) is a mm(G), 
which we write in short as mwm(G) = mm(G) and its weight is given 
by Wt(mwm(G)) = |mm(G)|. Similarly, if G is an undirected and weighted 
graph with Wt(e) = c for all edges e in G and c is a constant then also we 
have mwm(G) = mm(G) with the weight of the matching as Wt(mwm(G)) = 
cx |mm(G)|. Given a bipartite graph G, the problem of Mazimum Weight 
Bipartite Matching (MWBM) computes a MWM of G. 


1.2. Our Contribution 


In [19,20], Kao et al. proposed a decomposition theorem and algorithm 
for computing weight of a MWBM of the bipartite graph G. This paper 
revise the existing decomposition theorem and use it efficiently to design an 
improved version of the decomposition algorithm which estimates the weight 
of a MWBM of G in time O(,\/|V|W'/k(|V|, W’/N)) by taking algorithm 
designed by Feder and Motwani [11] as a base algorithm, where k(x, y) = 
log x/ log(x?/y). 

This algorithm bridges a gap between the best known time complexity 
of computing a Maximum Cardinality Matching (MCM) and that of com- 
puting a MWBM of a bipartite graph. In the best case, computation of 
weight of a MWBM takes O(./|V||E|/k(|V|,|£|)) time which is the same 
as the complexity of the Feder and Motwani’s algorithm [11] for comput- 
ing MCM of an unweighted bipartite graph; whereas in the worst case it 
takes O(\/|V|W/k(|V|, W/N)), that is, |E| < W’ < W. Further, we provide 
an interesting scaling property of the algorithm and a better bound of the 
parameter W’. However, it seems to be a challenging problem to get rid 
of W or N from the complexity. 

The modified algorithm works well for general W, but is best known 
for W’ = o(|E|log(|V|N)). We also design an algorithm to construct a 
minimum weight cover of a bipartite graph in O(./|V|W'/k(|V|, W’/N)) 
time. It identifies the edges involved in MWBM. It is also possible to 
use other algorithms as a subroutine, for example, algorithms by Hopcroft 
and Karp [17] and Alt et al. [1] in which case the running times of our 
algorithm will be O(,/|V|W’) and O((|V|/log|V|)!/2W’), respectively. An 
experimental evaluation of randomly generated bipartite graphs shows that 
the proposed improvement is significant in general. 
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1.3. Roadmap 


The rest of the paper is organized as follows. In Section 2, we give a detailed 
summary of existing maximum matching algorithms and their complexities 
for unweighted and weighted bipartite graphs. Section 3 describes a modified 
decomposition theorem and an algorithm to compute the weight of a MWBM. 
The complexity analysis of the algorithm is discussed in Section 4. The 
algorithm to compute a minimum weight cover of a bipartite graph is given 
in Section 5, which is used to find the edges of a MWBM. Section 6 provides 
the experimental comparisons between the modified algorithm and Kao et 
al.’s algorithm for randomly generated bipartite graphs. We summarize the 
results in Section 7. 


2 Survey of Maximum Matching Algorithms in 
Bipartite Graphs 


The problem of computing maximum matching in a given graph is one of 
the fundamental algorithmic problems that has played an important role in 
the development of combinatorial optimization and algorithmics. A survey 
of some of the well known existing maximum (cardinality) matching and 
maximum weight matching algorithms for the bipartite graph are summarized 
in Table 1 and Table 2, respectively. The algorithms with best asymptotic 
bound are indicated by “x” in these tables. A more detailed and technical 
discussion of the algorithms can be found in textbooks [21, 25, 26]. 


2.1 Maximum Cardinality Matching 


For unweighted bipartite graphs, Hopcroft-Karp [17] algorithm, which is 
based on augmenting path technique, offers the best-known performance 
for finding maximum matching in time O(|E|,/|V]). In the case of dense 
unweighted bipartite graphs, that is with |E| = O(|V|?), slightly better 
algorithms exist. An algorithm by Alt et al. [1] obtains a maximum 
matching in O(|V|!°,/|E|/log|V|) time. In the case of |E| = @(|V|?), 
this becomes O(|E|,/|V|/log|V|) and is also \/log|V|-factor faster than 
Hopcroft-Karp algorithm. This speed-up is obtained by an application of 
the fast adjacency matrix scanning technique of Cheriyan, Hagerup and 
Mehlhorn [3]. The algorithm proposed by Feder-Motwani [11] has the time 


complexity O(|E|,\/|V|/k(|V|,|#])), where k(x, y) = log x/log(x*/y). 
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Table 1: Complexity survey of maximum unweighted bipartite matching 
algorithms. 


Year Author(s) Complexity 
1973 * | Hopcroft and Karp [17] O(\|E|V/IV1) 

1991 Alt, Blum, Mehlhorn and Paul [1] | O(|V|'°./|E|/log |V]) 
1995 * | Feder and Motwani [11] O(|E|/IVI/k(IV], EI) 


2.2 Maximum Weight Bipartite Matching 


Several algorithms have also been proposed for computing maximum weight 
bipartite matching, improving both theoretical and practical running times. 
The well known Hungarian method, the first polynomial-time algorithm, 
was introduced by Kuhn [22] and Munkres [23]. Fredman and Tarjan [12] 
improved this with running time O(|V|(|E| +|V|log|V]|)) for a sparse graph 
by using Fibonacci heaps. An O(|V|?/4|E|log N)-time scaling algorithm 
was proposed by Gabow [13] under the assumption that edge weights are 
integers. A different and faster scaling algorithm was given by Gabow and 
Tarjan [14] with running time O(,/|V]||E|log(|V|N)). Kao et al. [20] pro- 
posed an O(\/|V|W/k(|V|, W/N))-time decomposition technique under the 
assumption that weights on the edges are positive and W = o(|E| log(|V|N)). 

In addition to the above exact algorithms, several randomized and 
approximate algorithms are also proposed, see for example [10,24]. For a 
tight lower bound for the weights of maximum weight matching in bipartite 
graphs, please refer to [6]. 


3 Refined Decomposition Theorem for MWBM 


We now propose a modified decomposition theorem which is a generalization 
of the existing decomposition theorem originally proposed by 
Kao et al. [19,20] and use it to develop a revised version of the decom- 
position algorithm to decrease the number of iterations and speed up the 
computation of the weight of a MWBM. Let G = (V = Vi UVa, E, Wt) be 
an undirected, weighted bipartite graph having V; and V2 as a partition of 
the vertex set V. Further, let EF = {e1,€2,...,e€z)} be set of edges of G 
with weights Wt(e;) = w; for 1 < i < |E|, where w1,we,..., w)z) are not 
necessarily distinct. As defined earlier, let N be the maximum edge weight, 
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Table 2: Complexity survey of maximum weight bipartite matching algo- 
rithms. 


Year(s) Author(s) Complexity 
1955, Kuhn [22], O(|V|*) 
1957 Munkres [23] (Hungarian method) 
1960 Tri [18, 25] O(\V|?|E]) 
1969 Dinic and Kronrod [9, 25] O(\V|) 
1984, 1987 « | Fredman and Tarjan [12] O(|V|(|EZ| + |V| log |V])) 
1985 Gabow [13] O(|V|3/4|B| log N) 
1989 * Gabow and Tarjan [14] O(r/|V]|E| log(|V|N)) 
1999 Kao, Lam, Sung and Ting [19] | O(,//V|W) 
2001 x Kao, Lam, Sung and Ting [20] | O(,/|V|W/k(|V|, W/N)) 
2014 x Das and Kapoor [7] O(/|V|W’) 
(This work) O((|V|/ log |V|)!/2W’) 

O(/IVIW'/k(\V|, W'/N)) 


that is, for all i € {1,2,...,|E]}, 0 < w; < N, and W = 2 w; be the 
1<t<|E| 
total weight of G. 

Our algorithm considers several intermediate graphs with lighter edge 
weights. During this process, it is possible that the weights of some of 
the edges may become zero. An edge e € E is said to be active if its 
weight Wt(e) > 0, otherwise, it is said to be inactive, that is when Wt(e) = 0. 
Let there be m’ (< |F|) distinct edge weights in the current working graph 
where wy < wa <-++: < Wm/_1 < Wm’. We denote the first two distinct 
maximum edge weights in the current working graph by H, and H2 (< Hj), 
respectively. Assign Hy = 0 in case m’ = 1. 


We first build two new graphs referred to as G;, and Ga from a given 
weighted bipartite graph G. For any integer h € [1, N], we decompose the 
graph G into two lighter weighted bipartite graphs G, and Ga as proposed 
by Kao et al. [19,20]. A minimum weight cover is a dual of maximum 
weight matching [20]. A cover of G is a function C: Vj U V2 > No such 
that C(v1) + C(ve) > Wie), Vor € Vi and Vu2g € Vo, where e = {v4, vo}. 
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Let Wi(C) = x C(x). A cover C is minimum weight cover if Wt(C) is 


rEV{UV2 
minimum. 


Formation of G, from G: The graph Gy, is formed by including those 
edges {u,v} of G whose weights Wt(u, v) lie in the range [VN —h+1, N]. 
Each edge {u,v} in graph Gy is assigned weight Wt(u,v) — (N — h). 
For illustration, G is constructed by the maximum weight edges of G 
and assigned unit weight to each edge. 


Formation of GA from G: Let Cy, be the minimum weight cover of Gp. 


The graph es is formed by including every edge {u,v} of G whose 
weight satisfies the condition 


Wi(u, v) — Ch(u) — Ch(v) > 0. 
The weight assigned to such an edge is Wt(u, v) — Ch(u) — Ch (v). 


Theorem 3.1 (The Decomposition Theorem [20]) Let G be an undi- 
rected, weighted bipartite graph. Then 


(a) for any integer h € [1, N], 


Wt(mwm(G)) = Wt(mwm(Gr)) + Wt(mum(GS)), 
(b) in particular (trivial), for h = 1, 
Wt(mwm(G)) = Wi(mm(G1)) + Wi(mwm(G4)). 


Note that Theorem 3.1(b) is derived from Theorem 3.1(a), since for h = 1, 
we have 


mwm(G1) = mm(G1) 


and 


Wt(mwm(G1)) = Wt(mm(G‘1)) = |mm(G4)|. 


Theorem 3.1(b) is used recursively in Algorithm 1 [20], to compute the 
weight of a maximum weight matching of the graph G. 


46 Shibsankar Das 


Algorithm 1 Kao et al.’s algorithm [20] to compute the weight of a MWBM. 


Input: A weighted, undirected, complete bipartite graph G with positive 
integer weights on the edges. 


Output: Weight of a maximum weight matching of G, that is, Wt(mwm(G)). 


Compute-MwmM(G) 
1: Construct G; from G. 
: Compute mm(G‘) and find a minimum weight cover C of G1. 
: Construct GS from G and C\. 
: if GS is empty, 
then return Wt(mm(G1)); 
: else return Wt(mm(G))+Compute-MwM(G). 


aa wn 


Remark 3.2 Observe that for arbitrary h € [1,N], mwm(G),) need not 
be equal to mm(G),), that is, we cannot always conclude that mwm(G),) = 
mm(Gp). 


One of our objectives is to investigate those values of h for which mwm(G),) 
is equal to mm(G),) apart from the trivial value of h as 1 in each iteration 
of Algorithm 1 to generate Gy having all its edge weights as 1. 

In order to get the speed up whenever possible, by decreasing the 
number of iterations whenever possible, we revise Theorem 3.1(b) and 
propose Theorem 3.3 which gives a domain of h € [1, N] where mwm(G),) = 
mm(G';,) and as a consequence of that we can write 


Wt(mum(Gp)) = Wt(mm(Gp,)) = h * |mm(Gp)]- 


It works for h = 1 and performs well especially when the largest edge weight 
differs by more than one from the second-largest edge weight in the current 
graph in a decomposition step during an iteration. 


Theorem 3.3 (The Modified Decomposition Theorem) The follow- 
ing equalities hold for any integer h € {1, H, — Hz] where Hy and Hz (< H;) 
are the first two distinct maximum edge weights of graph G, respectively. We 
assign Hz = 0 in case all edge weights are equal. 


(a) mwm(Gp) = mm(Gp), 
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(b) Wt(mwm(G)) = h* Wt(mm(Gp)) + Wt(mwm(GS)). 


Proof: The proof of the above statements are based on the construction 
of new graphs G), and GS from G and Theorem 3.1(a). 


(a) To prove that for any integer h where 1 < h < Hy — H2, mwm(Gh) = 
mm(G,) holds true, it is enough to prove the same for the maximum 
value? of h, that is, for h = H;—Hy. As specified earlier, the construction 
of Gp, is done by choosing those edges {u,v} of G that have weight 


Wi(u,v) € [N-—h+1,N] = [Ai — (Ai — Ae) +1, Ai] = [He +1, A]. 


Since H; € [Hz +1, Hi], Gp has only the heaviest edges of G and each 
such edge is assigned the same weight. Thus, mwm(Gp) = mm(Ghr) 
for h= Ay — Ad. 


(b) Observe that h € [1, H, — He] and [1, H, — He] C [1,N]. So, by us- 
ing Theorem 3.1(a) we have, VA € [1, Hi — Ho], 


Wt(mwm(G)) = Wt(mwm(G))) + Wi(mwm(G2)). 


Also by using Theorem 3.3(a), mwm(G),) = mm(G,) for all h € [1, Hi — 
Hp]. Weight of each edge* {u,v} in G;, is exactly Wt(u,v) — (N —h) = 
H, — (H; —h) =h. Therefore, 


Wt(mum(Gp)) = hx Wt(mm(G,)) = h * |mm(Ga)|- 
Hence, for any integer h € [1, H; — Hy], 


Wt(mwm(G)) = Wt(mwm(G),)) + Wi(mwm(G2)) 
= h* Wt(mm(G))) + Wimwm(GR)). 


This completes the proof. 


Remark 3.4 The equality mwm(Gp) = mm(Gp) in Theorem 3.3(a) is not 
true forh > H, — Hy andh<N. 


°For illustration, consider h = c, where 1 < c < H, — Ho. Then as per the formation 
of G, from G, G. is built by choosing those edges of G that have weight Wt(u,v) € 
[N — (c—1), N]. Since, c— 1 >0 and N € [N — (c— 1), N] for any c € [1, Hi — He], Ge 
has only the heaviest edges of G. For optimization, choose h = H, — He, the maximum 
possible value of h. 

4Only maximum weight edges of G are included in Gp. 


48 Shibsankar Das 


To show that for any h € [H; — Hz +1, N] the statement mwm(G,) = 
mm(G,) is not true, it is enough to show the same essentially for h = Hi — 
Hyo+1. Observe that h = H, —H2+1> 2, since Hy > Ho. According to the 
construction of Gp, it is formed by edges {u, v} of G whose weights Wt(u, v) € 
[N —h+ 1,N] = [Ay (A, Ay 4 1) t 1, Ay] = [H2, Hy], that is, G,, is 
built with the maximum weight edges and second maximum weight edges 
of G, because {H, Hz} € [H2, Hi]. The weight of each heaviest edge {u, v} 
of G in Gs, is exactly 


Wt(u,v) —-(N-—h) = Hy -(M -— hy =h 


which is greater than or equal to 2 and that of each second heaviest edge {u, v} 
of G in Ga is exactly 


Wt(u,v) — (N —h) = Ho - (A, —h) = (Ag —- A) +h=(1—-h)+h=1. 


Hence, mwm(G;,) 4 mm(G);,) for such a value of h. 


(a) (b) 
Figure 1: (a) An undirected bipartite graph G with positive integer weights 
on the edges. (b) Considering h = H; — Hz +1 = 6, Gj, is extracted, 
but mwm(Gp,) 4 mm(G)). 


Example 3.5 Consider the graph shown in Figure 1(a). Leth = H,—Hp2+1. 
So, h = H, —-H2+1=9-4+4+1=6. As shown in Figure 1(b), Gp is 
formed by the edges {u,v} whose weights Wt(u,v) € [N-—h+1,N] = 
[9 -6+ 1,9] = [4,9] and their respective calculated weights are 6 and 1. 
Hence, mwm(G),) 4 mm(Ga). 
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We use the modified decomposition Theorem 3.3 to design a recursive 
Algorithm 2 to compute the weight of a mwm(G). 


Algorithm 2 Compute weight of a maximum weight matching of G. 


Input: A weighted, undirected, complete bipartite graph G with positive 
integer weights on the edges. 


Output: Weight of a maximum weight matching of G, that is, Wt(mwm(G)). 


WtT-MwesM(G) 
1: Assume that initially Wt(mwm(G)) = 0. 
: Find h = H, — Ho from the current working graph G. 
: Construct G;, from G. 
: Compute mm(G),). 
Find minimum weight cover Cy, of Gp. 
: Construct Ga from G and Ch. 
: if GS is empty (that is, GS has no active edge) 
then return h * |mm(G,)|; 
: else return h « |mm(G))| + WT-MwBM(G9). 


CwoONIA THR ww 


Example 3.6 Consider the bipartite graph as shown in Figure 2(a). Algo- 
rithm 2 finds the weight of a MWBM in just two iterations, as the algorithm 
is designed for the best h in every invocation of WT-MWBM(), whereas 
algorithm by Kao et al. [20] requires 500 iterations because it considers h = 1 
in every invocation of Compute-MwM( ). 


The correctness of the algorithm follows from the construction of G;, 
and Ga and the modified decomposition Theorem 3.3. 


4 'The Complexity of the Modified Algorithm 


Let G=(V = V, UV, E, Wt) be the initial input graph and N denotes the 

maximum edge weight of G, that is, for all i € {1,2,...,|E|},O<wi< N 

and W = a w; is the total weight of G. Further, let {w ,..., wm} be 
1<i<|B| 

the set of distinct edge weights of G, where m’ < |E]. 
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(a) (b) (c) (d) 


Figure 2: (a) An undirected, weighted bipartite graph G with positive 
integer weights on the edges. In the current graph G, h = 495. (b) Gj, is 
extracted. (c) C), is the weighted cover of G;,. (d) G? is formed from G, 
and C;,. Compute WT-MwsM(G9). 


Based on the constructions of G;, and G, the modified decomposi- 
tion Theorem 3.3 and Algorithm 2, we can easily observe that in the worst 
case the maximum number of possible iterations of WT-MwBM() is N, 
when fh = 1 in each iteration in the current working graph. Whereas in the 
best case, all the edge weights of G are the same and so we will have h = N for 
the present decomposition. As a consequence, the algorithm will terminate 
in the first iteration itself. 

As the complexity analysis of Algorithm 2 is almost similar to that 
presented elsewhere [20], the details are available in Appendix 8 (see page 66). 
The algorithm takes O(,/|V|W’/k(|V|, W’/N)) time to compute the weight 
of a mwm(G) by using the algorithm by Feder and Motwani [11], as a 
subroutine. 

Let L; consists of edges of remaining G (after (i—1)-th iteration) whose 
weights reduce in Ga in i-th iteration. Also let there be p iterations, |; = |L;| 
fori =1,2,...,p < N and h; = Hj — HAjg in the i-th iteration, where Hj 
and Hig (< Hj) are the first two distinct maximum edge weights of the 
remaining graph G after the (7 — 1)-th iteration. 

From the detailed complexity analysis we have, I, hi +lgh2+---+lphp = 
W. Let ) +lo+---4+h = W’. Observe that, in worst case, if h; = 1 


P 
for all i € [1,p], then W’ = Sk = W. And in best case, if hy = N, 
i=1 


then W’ = |E|. Moreover, the parameter W’ is smaller than W, essentially 
when the largest edge weight differs by more than one from the second-largest 
edge weight in the current working graph in decomposition step during at 
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least one iteration of the algorithm. Therefore in the best case°, it requires 
O(\/|V||E|/k(|V|, |E|)) time and in worst case O(,/|V|W/k(|V|, W/N)), to 
compute the weight of a maximum weight matching. That is, |E| < W’ < W. 

This time complexity bridges a gap between the best known time 
complexity for computing a Maximum Cardinality Matching (MCM) of an 
unweighted bipartite graph and that of computing a MWBM of a weighted 
bipartite graph. In the best case, for computation of weight of a MWBM, 
Algorithm 2 takes O(,/|V]|E|/k(|V|,|E|)) time which is the same as the 
complexity of the Feder and Motwani’s algorithm [11] for computing MCM 
of the unweighted bipartite graph; whereas in the worst case it (Algorithm 2) 
takes O(,/|V|W/k(|V|, W/N)) time which is the same as the complexity of 
the Kao et al.’s algorithm [20]. However, it is very difficult and challenging 
to get rid of W or N from the complexity. This modified algorithm works 
well for general W, but is best known for W’ = o(|E| log(|V|N)). 


4.1 Some More Advantages: Scaling Up and Down, and 
GCD Properties 


Some other advantages of the modified decomposition algorithm are stated 
by the following propositions. Let G = (V, E, Wt) be an undirected, weighted 
bipartite graph, FE = {e1,€2,...,€)z)} be the set of positive integer weight 
edges with weights Wt(e;) = w; > 0 (where 1 < i < |E|), N be the 
maximum edge weight and W = > w; be the total weight of G. The 
1<i<|B| 

modified decomposition Algorithm 2 computes the weight of a maximum 
weight bipartite matching of G in O(,/|V|W'/k(|V|, W’/N)) time, where 
|E|<W'<W. 


Proposition 4.1 (Multiplicative Scaling Up Property) Let G be a 
new weighted bipartite graph constructed by multiplying a large constant 
a € N to each edge weight w; of the initial weighted bipartite graph G. 
Then for both the graphs G and CG the complexity of Algorithm 2 remains 
O(\/|V|W'/k(|V|, W'/N)) where |E| < W’ < W; whereas for the graph G, 
the complexity of Algorithm 1 becomes O(,/|V|aW/k(|V|, aW/N)). 


Proof: As mentioned in the detailed complexity analysis of Algorithm 2 
(described in Appendix 8, page 66), let L; consists of edges of the remaining 


°In the best case, all the edge weights of G are the same. So, the algorithm terminates 
in just one iteration and hence W’ = O(|E)). 
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graph G (left after (i — 1)-th iteration), whose weights reduce in GS in 
the i-th iteration of WT-Mwsm(). Assume that there be p iterations for 
Algorithm 2, 1; = |£;| for i= 1,2,...,p < N and h; = Hj, — Hig in the i-th 
iteration, where Hj; and Hjz (< Hj) are the first two distinct maximum 
edge weights of the remaining graph G after (i — 1)-th iteration. From the 
detailed complexity analysis we have, 


lyhy + lohg +--+ + 1yhyp = W and htbt-::-+L,=W'. 


Observe that for the new graph G , the number of iterations in Algorithm 
Wt-Mwem(G) still remains p and in the computation of WT-MwBM(G), 
L; consists of 1; number of edges of the remaining graph G (after (i — 1)-th 
iteration), whose weights reduce in GA in i-th iteration of WT-MwsM(). In 
this case, if hi, = H!, — H{, in the i-th iteration, where H/, and H}, (< Hj) 
are the first two distinct maximum edge weights of the remaining graph G 
after (i — 1)-th iteration, respectively, then 


/ / e 
hi = Hi, — Hip =a* Hi —a* Hyp =ah; where i=1,2,...,p, 


Iyhy + lahg +--+ + ph, = hah + aha +--+ + lpahp = aW, 


and 

h+tbt+-::-+L,=W'. 
Therefore, the modified Algorithm 2 will take O(,/|V|W’/k(|V|, W’/N)) 
time to compute the weight of a mwm(G) by using the algorithm by Feder 


and Motwani [11] as a subroutine; whereas time required for the Kao et. 
al.’s Algorithm 1 is O(\/|VjaW/k(|V|, aW/N)) time. 


That is, multiplication by an integer constant to all the weight of edges 
of a weighted bipartite graph G does not affect the time complexity of the 
modified decomposition algorithm for computing the weight of a MWBM of 
the bipartite graph. The following remark talks about a conditional scaling 
down property of the algorithm for the graph G. 


Remark 4.2 (Multiplicative Scaling Down Property) Let we scale 
1 

down each edge weights of G by multiplying a factor of — and get a new 
a 


graph ee where a is the Greatest Common Divisor (GCD) of the positive 
edge weights of G. Then the time complexity of Algorithm 2 for computing a 
MWEBM of both the graphs G and G remains same. 
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Though during the complexity calculation of Algorithm 2 we have stated 
a bound for W’ as: |E| < W’ < W, but the following proposition gives a 
much better bound of the parameter W’. 


Proposition 4.3 (GCD Property) Let G = (V, E, Wt) be an undirected, 
weighted bipartite graph and E = {e1,€2,...,€)p\} be the set of positive 
weight edges with weights Wt(e;) = w; > 0 for1<i<|E|. Further, let the 
GCD of the positive edge weights of G is denoted by GCD(wi, w2,... , WR); 
then 


W 
E\<w'< < W. 
| | > GCD(wi, we,.--, WRI) 
Proof: Without going into more detailed and repeated writing, as 


mentioned in the previous Proposition 4.1, we have 


Iyhy + lohg +--+ +1yhp = W, where h; = Hj — H; and 
htbt-:--+h,=W'. 


Let g = GCD(w1, w2,.--, W)z)). Observer that, in any iteration 7 (where 
i=1,2,...,p) both Hj, and Hj are divisible by g. Hence, according to the 
definition of hjs, each h; = Hi, — Hig is also divisible by the factor g. 


Therefore, W’=1,+lo+---+ ly 


< ly (hi /g) oP lo(h2/g) Sadek Ip(hp/g) 


WwW 
< (ljhy + lghg +--+ + [php)/g = a =. 


This completes the proof. 


4.2 Complexity Analysis by Considering 
Other Base Algorithms 


We also analyze the complexity of Algorithm 2 by considering the Hopcroft- 
Karp algorithm [17] and the Alt-Blum-Mehlhorn-Paul algorithm [1] as base 
algorithms. 


With Respect to the Hopcroft-Karp Algorithm: Hopcroft-Karp algo- 
rithm [17] presents the best known worst-case performance for getting a 
maximum matching in a bipartite graph with a runtime of O(,/|V||E)). 
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Hence the recurrence relation for running time of Algorithm 2 with 
respect to the Hopcroft-Karp algorithm is 


T(V|,W',N) = O(/|V |) + 7(V],W", N") 
and = T(|V|,0,0) =0 


Therefore, T(|V|,W’, N) 


= O(o/|V[h) + O(/V] lz) +--+ + O(\/V Tl) 
Pp 

=O (vii) = O(,/\V|W’). 
i=l 


With Respect to the Alt-Blum-Mehlhorn-Paul Algorithm: A_ bit 
better algorithm for a dense bipartite graph is the Alt-Blum-Mehlhorn- 
Paul algorithm [1] which is (log |V|)!/?-factor faster than Hopcroft- 
Karp algorithm for maximum bipartite matching. Hence the time 
complexity, with respect to the Alt-Blum-Mehlhorn-Paul algorithm as 
a base algorithm, is O((|V|/log|V|)!/2W’) and it is (log |V|)!/?-factor 
faster than the above case. 


5 Finding a Maximum Weight Matching 


Algorithm 2 computes only the weight of a mwm(G) of a given graph G. To 
find the edges of a mwm(G), we give a revised algorithm for constructing a 
Minimum Weight Cover (MWC) of G which is a dual of maximum weight 
matching. As mentioned before, a cover of G is a function C: Vj U V2 > No 
such that C(v1) + C(ve) > Wt(v1, v2) Vor € Vi and vg € Vo. Let Wt(C) = 

y C(x). We say C is minimum weight cover if Wt(C) is minimum. Let 


rEV{UV2 
C be a MWC of a graph G. 


Lemma 5.1 ([20]) Let C4 be any minimum weight cover of GS. If C is 
a function on V(G) such that for every u € V(G), C(u) = Cp(u) + CA(u), 
then C is minimum weight cover of G. 


Using this lemma we design an O(,/|V|W'/k(|V|, W’/N))-time revised 
algorithm to compute a MWC of G (Algorithm 3). The correctness of this 
algorithm is clear from Lemma 5.1 and the time complexity analysis is similar 
to that given in the previous section. 


A Modified Decomposition Algorithm for Maximum Weight Bipartite 
Matching and Its Experimental Evaluation 55 


Algorithm 3 Calculate a MWC C of G. 


Input: A weighted, undirected, complete bipartite graph G with positive 
integer weights on the edges. 


Output: A minimum weight cover C’ of G. 


Mwc(G) 

: Assume that initially Wt(mwm(G)) = 0. 

: Find h + Hy — Ho from the current working graph G. 
: Construct Gp, from G. 

Compute mm(G;,). 

Find minimum weight cover Cy, of Gp. 

: Construct Ge from G and Ch. 

: if GS is empty (that is, G has no active edge) 

then return C}; 

: else 

CA & Mwc(G4); 

return C, where C(u) = C;,(u) + CA(u) for all nodes u in G. 


aioe 
ee © 


Now as deduced by Kao et al. in [20], finding a maximum weight 
matching by using the given vertex cover takes O(,/|V||E|/k(|V|,|E])) time. 
Since |E| < W’ < W, so altogether O(,/|V|W’/k(|V|, W’/N)) time requires 
to find a MWBM of G. 


6 Experimental Evaluations 


Algorithm 2 is efficient because of the modified decomposition Theorem 3.3. 
In order to understand the practical importance of Algorithm 2 (in compari- 
son with Algorithm 1), we report experimental evaluations of the same for 
the pseudo-randomly generated weighted bipartite graphs. 


6.1 Implementation and Experimental Environments 


We have implemented both Algorithm 1 (Kao et al.’s algorithm [20]) and 
Algorithm 2 (our modified algorithm) in C++ and compiled them using g++ 
4.8.2-19ubuntul compiler. All the experiments have been performed on a 
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Desktop PC with an Intel® Xeon® (E5620 @ 2.40 GHz) Processor, 32.00 GB 
RAM and 1200 GB Hard Disk, running the Ubuntu 14.04.1 LTS (Trusty 
Tahr) 64-bit Operating System. 


6.2. Input Data Description and Its Randomness 


For a frame of a fixed number of vertices and fixed total weight, we have 
pseudo-randomly generated several distinct weighted bipartite graphs by 
assigning random (uniformly distributed) weight to the randomly (uniformly 
distributed) picked up edges of a bipartite graph G. The outputs of these 
experiments for an input bipartite graph G are: 


(a) the number of iterations of WT-MwsBmM() and Compute-MwmM() in 
Algorithm 2 and Algorithm 1, respectively, for the pseudo-random 
bipartite graph G and 


(b) total time taken by the respective algorithms to compute the weight of 
a MWBM of G. 


As mentioned in [5,8,16], the approximate parameterized string matching 
problem in Stringology under the Hamming distance error model is compu- 
tationally equivalent to the MWBM problem in graph theory. The input 
data relation between the above problems are: 


(a) length of the pattern is equal to the total weight of the bipartite 
graph, and 


(b) the alphabet size of the pattern is equal to the number of vertices 
in a subset of the partition of the vertex set of the corresponding 
bipartite graph. 


6.3. Experimental Results 


We have tested the respective algorithms with large pseudo-random input 
data sets. The details are given below. 


Experiment 6.1 The first experiment is done for a total of 260 pseudo- 
randomly generated bipartite graphs, each of its weight is fixed to 1000 units 
where the size of the vertex set of bipartite graphs varies from 2 to 52. 
Each numerical row of Table 3 is corresponding to 10 different pseudo- 
randomly generated bipartite graphs, each of whose size of the vertex set and 
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Table 3: Efficiency comparison between Algorithm 2 and Algorithm 1 for 260 
pseudo-randomly generated weighted bipartite graphs as considered in Exper- 


iment 6.1. For each row, an average of Iterations 4 


t standard deviation and 


an average of Time (Sec.) are reported for 10 pseudo-randomly generated 
bipartite graphs for both the algorithms. 


ie Neavtiees Weight | Algorithm 2 (Our Modified Algorithm) | Algorithm 1 (by Kao et al.) 
of Graph # Iterations Time (Sec.) # Iterations | Time (Sec.) 
2 000 1.00 + 0.09 0.000015 279.00 + 19.64 0.000514 
4 000 8.40 + 1.14 0.000050 422.60 + 23.43 0.002356 
6 000 8.20 + 1.30 0.000084 397.60 + 37.80 0.003769 
8 000 6.60 + 0.89 0.000104 402.204 43.21 0.005922 
10 000 7.20 + 1.09 0.000146 444.40+ 15.75 0.009566 
12 000 8.60 + 0.98 0.000214 408.60 + 39.11 0.010882 
14 000 7.00 + 0.84 0.00025 418.60 + 27.44 0.014451 
16 000 10.40 + 1.67 0.000405 455.20 + 15.43 0.019964 
18 000 7.20 + 0.97 0.000376 366.40+ 33.96 0.018605 
20 000 7.40 + 1.02 0.00041 400.20 + 19.38 0.024284 
22 000 8.60 + 1.15 0.000525 391.20 + 33.39 0.025592 
24 000 9.00 + 0.95 0.000670 391.20 + 55.27 0.031218 
26 000 8.40 + 1.35 0.000722 391.20 + 19.83 0.035947 
28 000 8.40+ 1.21 0.000811 391.20 + 31.30 0.039951 
30 000 10.204 0.99 0.001089 391.20+ 33.39 0.044228 
32 000 9.60 + 0.89 0.001156 391.20 + 39.99 0.048768 
34 000 9.40 + 0.94 0.001217 391.20 + 29.63 0.053661 
36 000 9.60 + 1.03 0.001411 391.20 + 26.90 0.058682 
38 000 11.00 + 1.42 0.001680 391.20 + 31.43 0.064456 
40 000 8.60 + 1.04 0.001492 391.20 + 43.25 0.070000 
42 000 11.80 + 1.32 0.002226 391.20+ 23.46 0.077403 
44 000 11.00 + 2.27 0.002284 391.20 + 56.37 0.083169 
46 000 13.40 + 1.96 0.002924 391.20 + 27.48 0.089702 
48 000 13.80 + 1.19 0.003254 391.20 + 19.83 0.097006 
50 000 13.00 + 2.23 0.003327 391.20+ 33.39 0.103550 
52 000 12.204 2.02 0.003342 391.20 + 26.90 0.110369 


weight of the graphs are fixed. Only for this experimental result, for each row, 


an average of Iterations 4 


t standard deviation and an average of Time (Sec.) 


are reported for 10 different pseudo-randomly generated bipartite graphs for 


both algorithms. 


For example, the numerical row corresponding to 4 Vertices’ equal 
to 30 reports the following. For the 10 distinct pseudo-randomly generated 
bipartite graphs, each of whose size of the vertex set is 30 and weight is 
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1000 units, an average number of iterations of WT-MwBM() and Compute- 
MwM() in Algorithm 2 and Algorithm 1 are 10.20 and 391.20, respectively; 
whereas average time taken by the respective algorithms to compute the weight 
of a MWBM are 0.001089 second and 0.044228 second. 


The next two experiments are done over the graphs corresponding to the 
pseudo-randomly generated strings over the DNA alphabet ©’ = {A,C,G,T} 
of different lengths. 


Experiment 6.2 In this experiment we have fired the size of the vertex set 
of each of the 62 pseudo-randomly generated bipartite graphs to 8 with 62 
different weights ranging from 10 units to 3050 units. Please refer to Table 4 
for more details. Unlike the previous experiment, each row reports the 
comparisons of the number of iterations and time required by Algorithm 2 
and Algorithm 1 to compute a MWBM of a pseudo-randomly generated 
bipartite graph with fixed size vertex and total weight. 


Table 4: Experimental result for the 62 pseudo-randomly generated weighted 
bipartite graphs as considered in Experiment 6.2. The number of vertices in 
the vertex set of each of the bipartite graphs is fixed to be 8, but the weight 
of the bipartite graphs varies from 10 units to 3050 units. 


ie Martieds Weight | Algorithm 2 (Our Modified Algorithm) | Algorithm 1 (by Kao et al.) 

of Graph # Iterations Time (Sec.) # Iterations | Time (Sec.) 
8 10 3.00 0.000121 5.00 0.000152 
8 50 4.00 0.000109 13.00 0.000229 
8 100 9.00 0.000260 38.00 0.000965 
8 150 6.00 0.000165 52.00 0.001257 
8 200 5.00 0.000141 81.00 0.001641 
8 250 5.00 0.000142 109.00 0.002703 
8 300 18.00 0.000419 133.00 0.003255 
8 350 5.00 0.000144 116.00 0.002648 
8 400 6.00 0.000192 139.00 0.003666 
8 450 4.00 0.000122 192.00 0.004288 
8 500 31.00 0.000745 189.00 0.004863 
8 550 6.00 0.000165 203.00 0.004828 
8 600 6.00 0.000159 277.00 0.006564 
8 650 6.00 0.000166 298.00 0.006679 
8 700 9.00 0.000195 186.00 0.003468 
8 750 8.00 0.000222 268.00 0.007218 
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Table 4 (cont.) 


4 Wentices Weight | Algorithm 2 (Our Modified Algorithm) | Algorithm 1 (by Kao et al.) 
of Graph # Iterations Time (Sec.) # Iterations | Time (Sec.) 
8 800 7.00 0.000188 243.00 0.005828 
8 850 8.00 0.000193 390.00 0.009581 
8 900 8.00 0.000205 380.00 0.010258 
8 950 11.00 0.000297 395.00 0.011187 
8 000 11.00 0.000248 368.00 0.009515 
8 050 12.00 0.000230 466.00 0.012499 
8 100 5.00 0.000135 535.00 0.012997 
8 150 6.00 0.000126 405.00 0.008078 
8 200 8.00 0.000205 397.00 0.008670 
8 250 6.00 0.000168 602.00 0.014516 
8 300 6.00 0.000168 648.00 0.016127 
8 350 13.00 0.000283 553.00 0.015214 
8 400 6.00 0.000158 639.00 0.016306 
8 450 10.00 0.000242 426.00 0.009439 
8 500 94.00 0.002554 456.00 0.011925 
8 550 9.00 0.000267 591.00 0.017001 
8 600 6.00 0.000162 775.00 0.016471 
8 650 7.00 0.000189 676.00 0.015674 
8 700 6.00 0.000162 665.00 0.013834 
8 750 6.00 0.000162 860.00 0.021617 
8 800 6.00 0.000163 701.00 0.017824 
8 850 6.00 0.000150 777.00 0.016954 
8 900 6.00 0.000142 827.00 0.019956 
8 950 8.00 0.000194 788.00 0.018098 
8 2000 6.00 0.000168 690.00 0.016825 
8 2050 6.00 0.000175 584.00 0.011463 
8 2100 37.00 0.000864 727.00 0.016040 
8 2150 5.00 0.000139 585.00 0.009586 
8 2200 8.00 0.000186 747.00 0.016389 
8 2250 6.00 0.000175 897.00 0.022177 
8 2300 10.00 0.000252 1091.00 0.028346 
8 2350 6.00 0.000164 740.00 0.014673 
8 2400 8.00 0.000185 954.00 0.022879 
8 2450 6.00 0.000168 998.00 0.023334 
8 2500 18.00 0.000439 735.00 0.017273 
8 2550 5.00 0.000159 1086.00 0.023379 
8 2600 7.00 0.000200 1035.00 0.027333 
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Table 4 (cont.) 


aervertioes Weight | Algorithm 2 (Our Modified Algorithm) | Algorithm 1 (by Kao et al.) 

of Graph # Iterations Time (Sec.) # Iterations | Time (Sec.) 
8 2650 6.00 0.000155 1313.00 0.032712 
8 2700 4.00 0.000133 1348.00 0.031344 
8 2750 10.00 0.000254 922.00 0.020268 
8 2800 9.00 0.000271 1030.00 0.028735 
8 2850 15.00 0.000361 1206.00 0.029624 
8 2900 42.00 0.000958 1144.00 0.026569 
8 2950 5.00 0.000153 1266.00 0.030598 
8 3000 6.00 0.000180 1249.00 0.033715 
8 3050 8.00 0.000216 1117.00 0.027538 


Experiment 6.3 In the final experiment also we have fixed the size of the 
vertex set of each of the 71 pseudo-randomly generated bipartite graphs to 
& with 71 different and large weights ranging from 1000 units to 700000 
units. See Table 5 for more details. Each row reports the comparison of 
the number of iterations and time taken by Algorithm 2 and Algorithm 1 to 
compute a MWBM of a pseudo-randomly generated bipartite graph with fixed 
size verter and weight. 


Table 5: Experimental result for the 71 pseudo-randomly generated bipartite 
graphs as considered in Experiment 6.3. The cardinality of the vertex set of 
each of the bipartite graphs is fixed to be 8, but the weight of the bipartite 
graphs varies largely from 1000 units to 700000 units. 


ae Wertieas Weight | Algorithm 2 (Our Modified Algorithm) | Algorithm 1 (by Kao et al.) 

of Graph # Iterations Time (Sec.) # Iterations | Time (Sec.) 
8 1000 11.00 0.000241 368.00 0.009156 
8 10000 7.00 0.000196 4284.00 0.009156 
8 20000 8.00 0.000193 7722.00 0.009156 
8 30000 6.00 0.000161 14942.00 0.352380 
8 40000 12.00 0.000326 19786.00 0.589580 
8 50000 8.00 0.000213 22172.00 0.619935 
8 60000 8.00 0.000207 28763.00 0.806167 
8 70000 8.00 0.000215 22042.00 0.588712 
8 80000 9.00 0.000190 28054.00 0.662195 
8 90000 10.00 0.000240 21440.00 0.447115 
8 100000 36.00 0.000970 32975.00 0.865868 
8 110000 10.00 0.000222 53322.00 1.320364 
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Table 5 (cont.) 


4 Nerticas Weight | Algorithm 2 (Our Modified Algorithm) | Algorithm 1 (by Kao et al.) 
of Graph # Iterations Time (Sec.) # Iterations | Time (Sec.) 
8 20000 8.00 0.000199 57741.00 466250 
8 30000 6.00 0.000156 54970.00 .209947 
8 40000 9.00 0.000204 50849.00 102937 
8 50000 12.00 0.000221 65220.00 502660 
8 60000 9.00 0.000183 44405.00 0.868532 
8 70000 8.00 0.000188 71264.00 -668321 
8 80000 9.00 0.000212 68565.00 391208 
8 90000 7.00 0.000182 58417.00 .261169 
8 200000 8.00 0.000199 84882.00 2.177765 
8 210000 21.00 0.000411 55160.00 0.861068 
8 220000 7.00 0.000171 90913.00 -931699 
8 230000 10.00 0.000247 87201.00 2.051072 
8 240000 15.00 0.000315 112402.00 2.794717 
8 250000 8.00 0.000180 102078.00 2.032557 
8 260000 17.00 0.000353 105322.00 2.745294 
8 270000 10.00 0.000251 88840.00 2.062840 
8 280000 15.00 0.000368 94300.00 2.191243 
8 290000 8.00 0.000191 79909.00 1.639328 
8 300000 18.00 0.000364 05128.00 2.443022 
8 310000 8.00 0.000188 20579.00 2.260459 
8 320000 8.00 0.000216 25597.00 3.134282 
8 330000 10.00 0.000236 51182.00 3.940926 
8 340000 9.00 0.000196 66492.00 3.507791 
8 350000 7.00 0.000181 51689.00 3.194340 
8 360000 7.00 0.000186 66928.00 3.850720 
8 370000 9.00 0.000212 67154.00 3.666068 
8 380000 10.00 0.000162 47368.00 2.885593 
8 390000 21.00 0.000388 37803.00 2.969532 
8 400000 9.00 0.000202 58026.00 3.690540 
8 410000 10.00 0.000256 53238.00 4.177567 
8 420000 15.00 0.000350 68440.00 4.543331 
8 430000 8.00 0.000174 95902.00 4.567199 
8 440000 7.00 0.000192 65922.00 4.078337 
8 450000 8.00 0.000199 83992.00 4.580142 
8 460000 8.00 0.000190 208746.00 4.302464 
8 470000 14.00 0.000234 229321.00 5.470295 
8 480000 8.00 0.000210 99475.00 5.379294 
8 490000 10.00 0.000266 239623.00 6.374744 
8 500000 11.00 0.000238 86026.00 4.691640 
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2b $eiHeea Weight | Algorithm 2 (Our Modified Algorithm) | Algorithm 1 (by Kao et al.) 
of Graph # Iterations Time (Sec.) # Iterations | Time (Sec.) 
8 510000 12.00 0.000298 201041.00 4.919859 
8 520000 9.00 0.000218 189501.00 4.317819 
8 530000 9.00 0.000205 151069.00 3.144941 
8 540000 8.00 0.000219 230280.00 5.777490 
8 550000 7.00 0.000164 254817.00 5.943550 
8 560000 20.00 0.000433 240510.00 6.015034 
8 570000 7.00 0.000178 260619.00 5.542112 
8 580000 7.00 0.000177 230100.00 4.817907 
8 590000 11.00 0.000239 240188.00 5.613709 
8 600000 8.00 0.000202 260924.00 5.703650 
8 610000 7.00 0.000179 261019.00 6.318097 
8 620000 13.00 0.000295 281111.00 5.626782 
8 630000 9.00 0.000227 276200.00 5.939048 
8 640000 10.00 0.000227 286193.00 7.329997 
8 650000 10.00 0.000213 255659.00 5.345079 
8 660000 6.00 0.000165 321451.00 8.061329 
8 670000 7.00 0.000170 243797.00 5.014862 
8 680000 8.00 0.000173 286228.00 6.904159 
8 690000 19.00 0.000396 306680.00 8.027268 
8 700000 4.00 0.000134 330990.00 7.337250 


7 Conclusions and Future Work 


We have fine-tuned the existing decomposition theorem originally proposed 
by Kao et al. in [20], in the context of maximum weight bipartite match- 
ing and applied it to design a revised version of the decomposition algo- 
rithm to compute the weight of a maximum weight bipartite matching in 
O(\/|V|W'/k(|V|, W'/N)) time by employing an algorithm designed by Feder 
and Motwani [11], as a base algorithm. We have also analyzed the algorithm 
by using the Hopcroft-Karp algorithm [17] and the Alt-Blum-Mehlhorn-Paul 
algorithm [1] as a base algorithm, respectively. 

The algorithm performs well especially when the largest edge weight 
differs by more than one from the second-largest edge weight in the cur- 
rent working graph during an invocation of WT-MwsM() in any itera- 
tion. Further, we have given a scaling property of the algorithm and a 


bound of the parameter W’ as |E| < W’ < W, 


=< 

~ GCD(wj, we,..- , W\B}) 
where GCD(w1, we, ..., W)z)) denotes the GCD of the positive edges weights 
{w1,W2,---,W)z)} of the weighted bipartite graph. The algorithm works well 
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for general W, but is the best known for W’ = o(|E|log(|V|V)). The exper- 
imental study shows that the performance of the modified decomposition 
algorithm is satisfactory. 

The modified algorithm seems much better on average, whereas its 
complexity is the same as that of the original algorithm in the worst case. 
Computing the average running time would certainly be an interesting 
future work. 
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8 Appendix: Detailed Complexity Analysis of 
Algorithm 2 


Here we give complexity analysis of Algorithm 2 in general. It is almost 
similar as done in the paper [20]. We assume that a maximum heap [4] is 
used to store the distinct edge weights along with the associated edges of G. 

Let the running time of WT-MwBM(G) be T(|V|, W’, NV) excluding the 
initialization. Let DL be the set of the heaviest weight edges in G. So up 
to the Step 3, construction of Gp, requires O(|L|log|E|) time. Step 4 takes 
O(\/|V|L|/k(|V],|L|)) time by using Feder and Motwani’s algorithm [11] 
to compute mm(G,). In Step 5, C;, can be found in O(|L|) time from this 
matching. Let L; be the set of edges of G adjacent to some node u with 
C,,(u) > 0, that is, Ly consists of edges of G whose weights reduce in Ge. 
Let 1, = |L|. Step 6 updates every edge of L, in the heap in O(I; log |E|) 
time. Since L C Ly, Step 1 to 6 takes O(4/|V|li/k(|V|, 11)) time altogether. 
Let 1; = |L;| fori = 1,2,...,p < N and h; = Hj — Hp for i-th phase of the 
recursion, where L; consists of edges of remaining G whose weights reduce 
in Ga on i-th iteration. Note that, 


Iyhy + lghe +---+1php = W. 
Let l) +l2 +--+: +l) = W’. Observe that if h; = 1 for alli € [1,p], 
P 
then W! = by = W. Step 7 uses at most T(|V|,W”, N”) time, where 
i=1 
W" (< W’) is the total weight of GS and N” (< N) is the maximum edge 
weight of fe Hence the recurrence relation for running time is 
T(|V|,W', N) = O(V|V|a/k(\V]4)) +TV|,W",N") and 
T(|V|,0,0) =0 
Therefore, T(|V|,W’, N) 


-0( 4) .0( vo) + | a - 
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_ A iy cae cl 
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2p we) >. Sie: 
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Let f(x) = xlogz. Note that it is a convex function, so by Jensen’s 
inequality®, 


ar 
. . i=1 
do Wilogli = 3° #4) 2 of | =— | = ef (“) 


/ / / Ww’ 
ere = Wie =0 (w"tog * ) 
p Pp p N 
This leads to the running time complexity as follows. 


T(|V|,W’,N) 


/ / 
-o( al (wrogivi?— wrt) 


log |V| 


Jiviw' 


2 
log |V|/log yprny 


= O(/|V|W'/k(|V], W'/N)). 


This is better than the O(,/|V|W/k(|V|, W/N)) time as mentioned in [20]. 

The parameter W’ is smaller than W which is the total weight of G, 
essentially when the heaviest edge weight differs by more than one unit 
from the second heaviest edge weight in a current working graph during 
a decomposition in any iteration of the algorithm. In the best case, the 
algorithm takes O(,/|V||E|/k(|V|, |E|)) time to compute a maximum weight 
matching and in worst case O(,/|V|W/k(|V|, W/N)), that is, |E| < W! < W. 
This time complexity bridges a gap between the best known time complexity 
for computing a Maximum Cardinality Matching (MCM) of an unweighted 
bipartite graph and that of computing a MWBM of a weighted bipartite 
graph. 

However, it is very difficult and challenging to get rid of W or N from 
the complexity. This modified algorithm works well for general W, but is 
best known for W’ = o(|E| log(|V|N)). 


°Jensen’s Inequality [15]. If f(x) is a convex function on an interval J and 


[41, H2,-+-, [mn are positive weights such that > pe = 1 then f (305 es) = Ye eee). 
i=0 
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